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ABSTRACT 


The  Unsteady  Reynolds- Averaged  Navier  Stokes  (URANS)  approach  applied 
in  conjunction  with  the  Eddy  Interaction  Model  (EIM)  is  used  here  to  con¬ 
duct  numerical  simulations  of  aerosol  particle  deposition  on  the  human  form. 
Numerical  simulation  of  a  two-phase  turbulent  impinging  jet  flow  is  studied 
to  validate  the  proposed  methodology.  The  numerical  deposition  results  are 
compared  with  some  relevant  experimental  measurements.  The  results  of 
this  validation  show  that  the  standard  EIM  with  turbulent  tracking  tends  to 
overpredict  the  deposition  efficiency.  Greatly  improved  results  were  achieved 
by  using  a  near- wall  correction  with  the  EIM.  After  the  validation  of  the 
methodology,  a  number  of  detailed  numerical  simulations  of  aerosol  particle 
deposition  on  a  human  form  were  conducted.  Simulations  were  conducted  for 
13  aerosol  particle  sizes  ranging  from  O.f  /im  to  40  //in.  For  each  of  these  par¬ 
ticle  sizes,  three  orientations  of  the  human  form  relative  to  the  incident  wind 
direction  were  simulated  (viz.,  with  the  human  form  oriented  0,  90  and  180 
degrees  with  respect  to  the  direction  of  the  incident  wind).  High-resolution 
isopleths  of  the  deposition  probability  on  the  surface  of  the  human  form  are 
presented.  To  provide  the  quantitative  information  for  future  applications  by 
researchers,  extensive  data  files  have  been  produced,  providing  information 
on  the  cell  area  for  each  triangle  used  in  the  tessellation  of  the  surface  of  the 
human  form,  the  associated  body  part  associated  with  this  triangular  cell, 
and  the  computed  deposition  probability  associated  the  triangular  cell. 
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Chapter  1 
Introduction 


This  study  is  concerned  with  the  prediction  of  the  deposition  of  aerosol  par¬ 
ticles  on  the  human  form.  The  methodology  utilizes  computational  fluid  dy¬ 
namics  (CFD)  simulations  to  determine  the  complex  turbulent  flow  around 
the  human  form,  and  this  information  is  used  as  the  input  to  a  Lagrangian 
stochastic  model  to  predict  the  motion  around  and  deposition  of  aerosol  par¬ 
ticles  on  the  human  form.  The  is  done  by  studying  the  statistics  of  a  large 
number  of  “heavy”  particles,  each  of  which  is  moved  stochastically  in  the 
turbulent  flow  field  (reconstructed  from  the  CFD  simulations)  around  the 
human  form. 

There  are  five  mechanisms  by  which  an  aerosol  particle  can  be  deposited 
on  a  surface;  namely, 

1.  Interception; 

2.  Inertial  impaction; 

3.  Brownian  diffusion; 

4.  Sedimentation  (or  gravitational  settling); 

5.  Electrostatic  attraction. 

Aerosol  particles  of  1-30  //in  in  size  have  inertia  (or  mass) ,  and  they  will  not 
follow  exactly  the  sharp  streamlines  in  the  vicinity  of  a  bluff  body.  There¬ 
fore,  deposition  can  occur  as  a  result  principally  of  interception  and  inertial 
impaction.  Brownian  diffusion  is  only  important  if  the  particle  size  is  much 
less  than  1  pm.  Owing  to  the  fact  that  only  non-volatile  aerosol  particles 
are  considered  herein,  the  gravitational  settling  effect  is  included  implicitly 
in  our  Lagrangian  stochastic  model  (which  will  be  described  in  Section  2). 
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The  electrostatic  attraction  effect  on  particle  deposition  will  be  ignored  in 
this  study. 

There  are  three  basic  steps  which  are  required  to  study  aerosol  deposition 
on  the  human  form: 

1.  construction  of  an  accurate  grid  (mesh)  to  represent  the  human  form; 

2.  computation  of  the  complex  flow  fields  around  the  human  form  us¬ 
ing  the  Reynolds-Averaged  Navier  Stokes  (RANS)  or  unsteady  RANS 
(URANS)  approach  with  various  inflow  conditions; 

3.  development  and  implementation  of  the  Lagrangian  stochastic  model 
for  the  prediction  of  the  deposition  of  aerosol  particles  on  the  human 
form  for  various  cases,  and  collection  of  the  information  to  determine 
the  deposition  density  pattern  on  the  surface  of  the  body. 

The  remainder  of  this  report  is  structured  as  follows.  The  methodology 
to  achieve  the  three  basic  steps  above  is  discussed  in  detail  in  Section  2. 
In  Section  3,  three  test  cases  are  used  to  validate  the  proposed  numerical 
methodology.  In  Section  4,  the  results  for  the  detailed  numerical  simulations 
for  aerosol  particle  deposition  on  the  human  form  are  reported.  Finally, 
conclusions  are  given  in  Section  5. 
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Chapter  2 

Numerical  methods 


2.1  Numerical  methods  for  fluid  flow 


The  time- averaged  continuity  equation  and  Reynolds- Averaged  Navier  Stokes 
(RANS)  equations  for  an  incompressible  flow  are  given  as  follows: 


duj 

dxi 


0, 


(2.1) 


dui  dujUi  dp  d2Ui  d  (~r~r 

dt  dxj  dxi  dx 2  dxj  \ 


(2.2) 


where  the  Reynolds  averaging  of  a  quantity  is  denoted  using  an  overbar. 
Summation  is  implied  by  the  repeated  indices.  Here,  ut  and  u%  are  the  mean 
and  fluctuating  (turbulent)  velocities  in  the  ^-direction  (?'  =  1,  2,  3),  respec¬ 
tively. 

Simulations  using  the  RANS  equations  require  less  memory  and  compu¬ 
tational  power  than  other  numerical  methods  such  as  direct  numerical  sim¬ 
ulation  (DNS)  and  large-eddy  simulation  (LES).  However,  all  scales  of  the 
turbulence  field  need  to  be  modeled  (rather  than  simulated)  in  the  RANS  ap¬ 
proach.  Moreover,  the  averaging  procedure  in  the  RANS  equations  results  in 
additional  unknown  quantities  containing  products  of  fluctuating  velocities 
known  as  Reynolds  stresses  given  as  —upiy  In  order  to  parametrize  these 
Reynolds  stresses,  a  turbulence  closure  model  is  required.  The  Reynolds 
stresses  can  be  determined  by  a  turbulence  model,  either  through  the  turbu¬ 
lent  viscosity  hypothesis  (i.e.,  eddy  viscosity  approximation)  or  more  directly 
from  the  modeled  Reynolds  stress  transport  equations.  In  this  study,  the 
standard  k-e  model  [21]  is  used  to  parametrize  the  various  Reynolds  stresses. 

The  STREAM  code  [14]  is  used  in  this  study  to  solve  the  governing  equa¬ 
tions  (2.1)  and  (2.2)  to  generate  the  numerical  fluid  flow  fields  about  the  hu¬ 
man  form.  STREAM  is  a  fully-conservative,  block-structured  finite-volume 
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code  for  computational  fluid  dynamics,  which  employs  a  fully-collocated  stor¬ 
age  arrangement  for  all  transported  properties,  including  various  turbulence 
quantities  (e.g.,  turbulence  kinetic  energy,  viscous  dissipation  rate,  etc.). 
Within  an  arbitrary  non-orthogonal  coordinate  system,  the  velocity  vector 
is  decomposed  into  its  Cartesian  components,  and  these  are  the  components 
to  which  the  momentum  equations  relate.  Advective  cell-face  fluxes  are  ap¬ 
proximated  by  the  Upstream  Monotonic  Interpolation  for  Scalar  Transport 
(UMIST)  scheme  [15],  obtained  by  formally  imposing  Total  Variation  Di¬ 
minishing  (TVD)  constraints  on  Leonard’s  third-order  accurate  Quadratic 
Upstream  Interpolation  for  Convective  Kinematics  (QUICK)  scheme  [16]. 
A  second-order  fully  implicit  three-level  scheme  is  used  to  treat  the  tran¬ 
sient  (or,  local  tendency)  term.  The  mass  continuity  is  enforced  by  solving 
a  pressure-correction  equation  using  the  Semi-Implicit  Method  for  Pressure- 
Linked  Equations  (SIMPLE)  algorithm,  which  steers,  as  part  of  the  iterative 
sequence,  the  pressure  towards  a  state  in  which  the  mass  residuals  in  all  cells 
of  the  flow  domain  are  negligibly  small. 

All  transport  equations,  including  the  mean  momentum  and  turbulence 
equations,  are  discretized  and  solved  sequentially  as  part  of  the  SIMPLE 
algorithm  [17].  The  linearized  and  discretized  equations  obtained  in  the 
outer  iterations  of  the  SIMPLE  algorithm  are  solved  numerically  using  very 
efficient  iterative  linear  equation  solvers,  such  as  Stone’s  Strongly  Implicit 
Procedure  (SIP3D)  [18]  or  the  Conjugate  Gradient  Stabilized  (CGSTAB) 
method.  In  conjunction  with  a  fully-collocated  approach,  the  SIMPLE  al¬ 
gorithm  is  known  to  provoke  checkerboard  oscillations  in  the  pressure  field, 
reflecting  a  state  of  velocity- pressure  decoupling.  To  avoid  this,  the  widely 
used  method  of  Rhie  and  Chow  [19]  has  been  adopted  to  interpolate  the  cell- 
face  velocities  from  the  adjacent  nodal  values.  This  interpolation  essentially 
introduces  a  fourth-order  “pressure  smoothing”  to  remove  the  checkerboard 
oscillations  in  the  pressure  field.  Physical  diffusion  fluxes  are  approximated 
using  a  conventional  second-order  accurate  central  differencing  approach. 

For  mesh  generation,  the  “ray-casting”  approach  [20]  is  used  to  determine 
whether  a  computational  cell  lies  inside  the  complex  geometry  of  objects 
encountered  in  the  current  problem,  such  as  the  human  form  placed  within 
the  computational  domain.  If  the  cell  centroid  is  inside  an  object,  the  flag 
associated  with  this  cell  will  be  set  to  OBJECT.  Otherwise,  the  cell  flag  is 
set  to  FLUID,  allowing  an  efficient  matrix  solver  (e.g.,  SIP3D)  to  be  utilized 
for  the  solution  of  the  discretized  equations. 
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2.2  Particle  modeling  in  RANS  approach 

Although  DNS  and  LES  have  been  used  to  study  particle  deposition  in  simple 
turbulent  channel  or  pipe  flows  [1,  2,  7],  they  are  still  too  time-consuming 
to  be  used  for  complex  engineering  problems  involving  particle  deposition 
on  complex  surfaces.  As  an  alternative,  the  unsteady  Reynolds-Averaged 
Navier-Stokes  (URANS)  approach  is  used  here  for  the  prediction  of  the 
ensemble-averaged  flow  field,  which  is  used  in  conjunction  with  eddy  inter¬ 
action  models  (EIMs,  or  random  walk  models)  for  the  reconstruction  of  the 
instantaneous  flow  field.  EIM  is  based  on  a  Lagrangian  stochastic  approach 
where  individual  particles  are  allowed  to  interact  successively  with  discrete 
eddies.  Each  eddy  has  length,  velocity  and  lifetime  characteristic  scales  ob¬ 
tained  from  the  ensemble-averaged  flow  and  turbulence  fields  predicted  by 
the  RANS  approach.  The  end  of  the  interaction  between  a  particle  and  an 
eddy  occurs  when  the  time  of  interaction  between  the  particle  and  eddy  ex¬ 
ceeds  the  lifetime  of  the  eddy  or  when  the  particle  trajectory  crosses  the  eddy. 
At  this  instant,  a  new  interaction  of  the  particle  with  another  eddy  begins. 
The  trajectory  of  this  particular  particle  is  obtained  by  solving  the  equations 
of  motion  for  the  particle  [4],  In  simple  EIMs,  the  local  fluctuating  velocities 
are  obtained  by  multiplying  the  root-mean-square  fluid  fluctuating  velocity 
by  random  numbers  generated  from  a  Gaussian  probability  density  function 
with  zero  mean  and  unit  standard  deviation  at  the  start  of  an  eddy-particle 
interaction. 

Due  to  its  simplicity,  RANS(URANS)/EIM  is  widely  used  in  many  en¬ 
gineering  applications,  such  as  numerical  modeling  of  aerosol  deposition  in 
the  mouth-throat  region  in  human  upper  respiratory  system  [5].  However, 
limitations  associated  with  EIM  include: 

1.  particles  with  size  comparable  to  an  eddy  size  cannot  be  accurately 
treated; 

2.  the  motion  of  a  particle  does  not  account  for  the  effects  of  neighboring 
particles; 

3.  only  “one-way  coupling”  with  the  surrounding  flow  field  is  considered. 

Despite  such  deficiencies,  EIM  is  still  an  attractive  tool  due  to  its  compu¬ 
tational  efficiency  compared  to  DNS  and  LES,  and  is  especially  suitable  for 
complex  geometries.  Therefore,  URANS/EIM  is  selected  in  this  study  as 
the  principal  numerical  methodology  for  predicting  aerosol  dispersion  and 
deposition  on  the  complex  geometry  associated  with  the  human  form  under 
turbulent  flow  conditions.  Details  of  EIM  can  be  found  in  [6,  7],  and  its 
mathematical  formulation  is  summarized  briefly  below. 
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The  particle  equation  of  motion  for  a  sphere  can  be  written  as 


(g*^)  Pp  =  Fs  +  Fj  +  Fb  +  F„  +  F„.  (2.3) 

Here,  F(/  is  the  gravitational  acceleration  force  given  by 

Fs  =  (pP-p/)g;  (2-4) 

Fd  is  the  drag  force  which  we  parametrize  using  Stokes’  linear  resistance 
formulation  as 


F^  is  the  Basset  force  given  by 

Fb  =  67 \prp  f  K  (t  -t) 
Jtn 


cLt 


Up); 

(2.5) 

—dr; 

(2.6) 

F„  is  the  added  mass  or  the  inertia  of  the  virtual  mass  force  modeled  as 


2  3  D(uf-ur) 

”  3  pP/  Dt 


and,  Fa  the  force  caused  by  the  fluid  acceleration  given  by 

4  o  Du# 

F“ = rr-Pr ct 


(2.7) 


(2.8) 


In  Eqs.  (2.3)-(2.8),  rp  is  the  radius  of  the  spherical  aerosol  particle,  p  is  the 
particle  density,  g  is  the  gravitational  acceleration,  u  is  the  (instantaneous) 
velocity,  p,  is  the  fluid  viscosity,  K  (t  —  r)  is  a  kernel  function,  is  the 

material  derivative,  and  Rcp  =  P/IU/  ”pK2rp)  js  |-|ie  particle  Reynolds  number. 
The  subscripts  /  and  p  denote  the  fluid  and  particle  phases,  respectively. 
Except  for  the  drag  and  the  gravitational  acceleration  forces,  all  other  terms 
in  Eq.  (2.3)  are  only  significant  when  the  fluid  density  is  comparable  to,  or 
greater  than,  the  particle  density  [8].  In  all  cases  to  be  investigated  in  this 
report,  Pf/pp  <§C  1  is  assumed. 

By  reference  to  Figure  2.1,  let  us  define  a  turbulent  eddy  with  length 
scale  Le,  velocity  scale  ue  and  time  scale  re  as  follows  (here,  the  subscript 
“e”  denotes  “eddy”): 


ue  =  \l  -k,  based  on  an  isotropy  assumption, 


(2.9) 
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Particle 


Figure  2.1:  Basic  concept  underpinning  an  eddy  interaction  model  (EIM). 


Le  =  (C,f^—  (2.10) 

where  Le  is  the  radius  of  a  turbulent  eddy  depicted  by  a  circle  in  Figure  2.1. 
From  Eqs.  (2.9)  and  (2.10),  we  can  derive  the  eddy  time  scale  (or  lifetime  of 
a  turbulent  eddy)  as 


Te  = 


Le 

ZLe 


(C„)3/4  (k^/E) 


(2.11) 


Note  that  k  and  e  in  Eqs.  (2.9)-(2.11)  are  obtained  from  the  standard  k-s 
model,  and  C)(  =  0.09  is  a  turbulence  closure  constant  that  defines  the  eddy 
viscosity: 

IM  =  cj^.  (2.12) 

In  EIM,  all  particles  are  independent  of  each  other  and  there  is  no  particle 
collision;  viz.,  the  classic  dilute  gas  approximation  is  assumed.  Any  particle 
released  from  the  inlet  into  the  flow  field  is  allowed  to  interact  successively 
with  various  eddies,  each  with  a  radius  of  Le  obtained  from  the  pre-computed 
turbulence  (viz.,  k  and  s)  fields.  This  approach  was  referred  to  as  “one-way 
coupling”  earlier.  The  end  of  the  interaction  between  a  particle  and  an 
eddy  occurs  when  the  interaction  time  tn  exceeds  the  lifetime  of  an  eddy 
(viz.,  ti  >  re),  or  when  the  particle  crosses  the  boundary  of  the  eddy  (i.e. , 
when  the  particle  is  outside  of  the  eddy  depicted  by  a  circle  with  a  radius 
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of  Le  in  Figure  2.1).  At  this  instant,  a  new  interaction  for  the  particle 
with  another  eddy  is  started,  and  the  interaction  time  tt  is  reset  to  zero  to 
signal  the  beginning  of  another  interaction  of  the  particle  with  a  new  eddy. 
Furthermore,  Le  and  re  for  the  new  eddy  is  re-calculated  based  on  the  k-s 
model  solutions  at  the  particle  location  corresponding  to  the  event  U  =  0. 
The  particle  will  start  to  follow  another  trajectory  through  this  new  eddy  in 
accordance  to  its  equations  of  motion  as  described  by  Eq.  (2.3)  [4,  9]. 

As  implied  in  Eqs.  (2.3)-(2.8),  the  particle  is  driven  by  different  forces, 
which  are  related  to  the  instantaneous  fluid  velocity  vector  u/  =  u  =  rq. 
However,  the  RANS  approach  only  provides  the  mean  value  for  this  fluid 
velocity,  namely  11/  =  u  =  u*.  Therefore,  the  instantaneous  velocity  will 
need  to  be  “reconstructed”  based  on  the  mean  velocity  and  the  turbulence 
kinetic  energy,  k,  as  follows: 

Ui  =  ui  3  <q.  *  =  1,2,3,  (2-13) 

[(«!,  u2,  Us)  =  (u,  v,  w)]  where  the  components  of  the  fluctuating  velocity  are 
determined  from 


u 


v 


w 


(2.14) 


where  (Nu,  Nv,  Nw)  are  random  variates  drawn  from  a  Gaussian  probability 
density  function  with  zero  mean  and  unit  standard  deviation.  In  order  to 
obtain  deposition  density  statistics,  hundreds  of  thousands  to  millions  of 
particles  must  be  released  into  the  flow. 

Let  us  re-write  Eq.  (2.3)  ignoring  the  Basset,  virtual  mass  and  fluid  ac¬ 
celeration  forces  (owing  to  the  fact  that  p//pp  <C  1  :  for  the  cases  of  interest 
in  this  study): 


dup 

dt 


=  1 


P£ 

Pp 


g  + 


\PP3 


nrS, 


<f>( u  -  up) 


dup 

dt 


l -P-l 

Pp 


g+  ~ 


Tr 


u  - 


(2.15) 


where  0  =  (l  +  0.15  Re°'687)  [9]  and  rp  =  (ppdp)/(l8/i).  Eq.  (2.15)  when 
written  in  component  form  then  becomes 


dup 

dt 


Up) ) 


(2.16) 
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(2.17) 


dvp 

dt 


=  (I) c"  “  ’ 


-for  /  1  _  P/  ,  +  ,  x  ,  (ro  _  ^ 


dt 


Pp 


Tr 


(2.18) 


The  particle  equations  of  motion  describing  each  particle-eddy  interaction 
can  be  solved  numerically  by  using  a  first-order  explicit  time  stepping  scheme: 


dup 

dt 


u; 


71+1 


At 


_+  =  /  i  _  Pf 


Pp 


g  +  —  (u  -  uP) 


find  u”+1; 


(2.19) 


x7l+l  _  x7l 

p  —  p-  =  u”+1  find  x™+1.  (2.20) 

This  numerical  scheme  is  applied  to  determine  the  trajectory  for  each  parti¬ 
cle  released  into  the  flow.  The  particle  interacts  with  a  sequence  of  random 
eddies  representing  the  instantaneous  flow  (or,  velocity)  field,  and  the  trajec¬ 
tory  of  this  particle  is  tracked  until  it  collides  with  a  solid  surface  (e.g.,  the 
ground  surface,  the  surface  representing  the  human  body,  etc.)  or  exits  from 
the  boundaries  of  the  computational  domain.  The  time  step  that  is  used 
in  Eqs.  (2.19)-(2.20)  is  an  adaptive  time  step  and  reset  at  the  start  of  each 
eddy-particle  interaction  as  the  particle  enters  a  new  eddy.  The  time  step  is 
chosen  to  be  the  minimum  of  one  fifth  of  the  eddy  lifetime  or  one  fifth  of  the 
minimum  side  length  of  the  local  element,  Lmin,  divided  by  the  particle  local 
velocity  [10];  viz., 


At  =  min 


(  1  \  -^min 

w  + 55 + 


(2.21) 


Upon  a  particle-solid  surface  impaction  (i.e.,  a  particles  is  assumed  to 
deposit  and  stick  to  a  solid  surface  when  the  centroid  location  of  the  particle 
reaches  a  distance  from  the  surface  equal  to  the  radius  of  the  particle  or 
less),  the  particle  position  is  recorded  and  the  particle  is  assumed  to  be 
permanently  removed  from  the  flow.  In  other  words,  a  “sticky  wall”  condition 
is  imposed  so  that  no  rebounding  or  resuspension  of  a  particle  from  a  surface 
is  considered  in  the  modeling.  By  counting  the  number  of  particles  entering 
and  leaving  the  sample  section,  the  deposition  density  distribution  can  be 
obtained.  Near  a  solid  surface,  the  eddy  lifetime  is  assumed  to  be  a  constant 
and  equal  to  the  minimum  Lagrangian  time  scale  [11].  This  approach  is 
efficient  in  circumventing  the  problem  that  occurs  near  walls  (and  other  solid 
surfaces),  where  the  eddy  lifetimes  could  tend  to  vanishingly  small  values, 
which  will  increase  significantly  the  overall  computational  time. 
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Chapter  3 

Numerical  validation 


In  order  to  validate  the  implementation  of  EIM  and  assess  the  performance 
of  using  the  k-s  model  for  provision  of  the  flow  held  to  drive  the  dispersion  of 
the  aerosol  particles,  we  consider  a  two-phase  turbulent  impinging  jet,  which 
was  investigated  experimentally  in  [12]  and  numerically  in  [10].  A  schematic 
of  the  computational  domain  from  [12]  is  given  in  Figure  3.1  for  illustrative 
purposes. 

To  be  consistent  with  the  experiment  of  [12],  incompressible  and  steady 
state  fluid  how  at  a  standard  temperature  and  pressure  (STP)  will  be  as¬ 
sumed  for  the  simulations.  The  inlet  conditions  consist  of  a  constant  inlet 
velocity  Uj  =  10.5  m  s-1,  a  turbulence  intensity  of  5%  of  the  mean  velocity 
and  a  turbulence  length  scale  of  10%  of  the  inlet  diameter.  The  Reynolds 
number  based  on  the  inlet  diameter  is 


Re  = 


pUjD 

b 


(1.18)  (10.5)  (0.015) 
1.824  x  lO"5 


=  10189  ~  104. 


(3.1) 


The  standard  no  slip  boundary  condition  (BC)  will  be  applied  to  the  walls 
and  a  zero  gauge  pressure  BC  or  a  zero  gradient  (Neumann)  BC  combined 
with  bulk  mass  how  correction  will  be  used  at  the  outlet  of  the  computational 
domain. 

The  particle  Stokes  number  is  dehned  as 


Stk  =  C, 


PpUjdp 

9  pD 


(1) 


(1050) (10.5) (5  x  10“6)2 
(9)  (1.824  x  10“5)  (0.015) 


=  0.112, 


(3.2) 


where  Cc  ~  1  is  the  Cunningham  slip  factor,  pp  =  1050  kg  m-3  is  the 
particle  density,  and  dp  =  5  /.an  is  the  particle  diameter.  Experimental  data 
are  available  for  L/D  =  2,  4  and  6  9  (see  Figure  3.1). 

Because  the  computational  domain  is  axi-symmetric,  we  use  a  two-dimen¬ 
sional  axi-symmetric  code  to  solve  the  governing  equations  for  the  huid  how. 
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Figure  3.1:  A  schematic  of  the  computational  domain  for  the  two-phase 
impinging  jet  based  on  the  experimental  study  conducted  in  [12]. 


A  total  of  10, 000  particles  were  released  into  the  computational  domain. 
The  initial  velocity  of  each  particle  was  set  equal  to  the  fluid  velocity.  The 
particles  were  randomly  and  uniformly  distributed  in  space  at  the  inlet.  Two 
types  of  particle  tracking  were  tested:  namely,  “turbulent  tracking”  and 
“mean  flow  tracking”.  In  turbulent  tracking,  the  instantaneous  velocity  of 
the  flow  in  the  particle  equations  of  motion  is  taken  to  be  the  sum  of  the  mean 
and  fluctuating  flow  velocities.  In  mean  flow  tracking,  the  instantaneous 
velocity  of  the  flow  is  simply  taken  to  be  equal  to  the  mean  flow  velocity 
(viz.,  the  fluctuating  flow  velocity  is  assumed  to  be  zero). 

The  particle  deposition  density  profiles,  using  turbulent  tracking  (i.e., 
turbulent  dispersion)  and  mean  flow  tracking  (i.e.,  without  turbulent  disper¬ 
sion)  versus  the  radial  distance  from  the  stagnation  point,  r/D ,  are  shown 
in  Figures  3.2,  3.3  and  3.4  for  L/D  =  2,  4  and  6,  respectively.  These  particle 
simulations  were  performed  using  EIM.  The  experimental  data  of  Burwash 
et  al.  [12]  for  L/D  =  2,  4  and  6  are  shown  for  comparison.  Table  3.1  shows 
the  deposition  efficiency  defined  as  the  percentage  of  the  number  of  particles 
deposited  on  the  wall  for  the  three  cases.  The  experimental  results  and  nu¬ 
merical  results  from  CFX  5.7.1  [12]  are  summarized  in  the  table  as  well.  A 
perusal  of  Figures  3.2,  3.3  and  3.4  and  of  Table  3.1  shows  that  the  turbulent 
tracking  tends  to  overpredict  the  deposition  density  when  compared  to  the 
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L/D=2 


Figure  3.2:  Deposition  density  obtained  with  turbulent  tracking  and  mean 
flow  tracking  for  L/D  =  2  using  the  classical  EIM. 


experimental  measurements  for  L/D  =  2,  4  and  6.  This  overprediction  is 
not  unexpected  because  of  the  isotropy  assumption  employed  in  the  k-e  tur¬ 
bulence  model.  The  isotropy  assumption  gives  larger  turbulence  intensities 
(i.e.,  excessive  turbulence  kinetic  energy)  than  the  actual  component  normal 
to  the  wall,  which  plays  an  important  role  for  particle  deposition,  and  thus 
causes  more  deposition.  From  the  figures,  it  can  be  deduced  that  the  accu¬ 
rate  modeling  of  the  near- wall  effects  on  turbulence  features  is  important  for 
the  prediction  of  the  deposition  rates. 

In  contrast  to  the  turbulent  tracking,  the  mean  flow  tracking,  where  the 
particles  are  deposited  by  the  means  of  the  inertial  impaction  mechanism, 
underpredicts  the  deposition  density,  when  compared  to  the  experimental 
measurements.  The  underprediction  (for  L/D  =  2,  4  and  6)  of  the  deposition 
density  is  expected  because  at  a  size  of  5  /im,  most  of  the  particles  owing 
to  their  small  inertia,  will  tend  to  follow  the  mean  flow  field  and  only  a 
few  of  them  will  deposit  on  the  impingement  wall.  These  results  therefore 
suggest  that  the  turbulent  dispersion  mechanism  which  drives  the  turbulent 
tracking  is  the  main  mechanism  that  causes  the  particles  to  deposit  on  the 
impingement  wall. 

In  a  turbulent  flow,  the  turbulent  dispersion  (i.e.,  by  the  fluid  flow  fluc¬ 
tuations)  is  the  main  mechanism  for  particle  deposition,  especially  for  small 


21 


UD=4 


Figure  3.3:  Deposition  density  obtained  with  turbulent  tracking  and  mean 
flow  tracking  for  L/D  =  4  using  the  classical  EIM. 


L/D=6 


Figure  3.4:  Deposition  density  obtained  with  turbulent  tracking  and  mean 
flow  tracking  for  L/D  =  6  using  classical  EIM. 
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L/D 

Turbulent  tracking 

Mean  flow  tracking 

CFX  5.7.1 

Exp.  Burwash  (2006) 

2 

56% 

1.4% 

75% 

16.5% 

4 

57% 

2.9% 

77.4% 

17% 

6 

41% 

2.3% 

82% 

17.8% 

Table  3.1:  Comparison  of  the  deposition  efficiency  obtained  from  turbulent 
tracking  and  mean  flow  tracking  with  experimental  measurements. 


particle  sizes  (5  /mi  and  Stk  =  0.11).  Therefore,  the  appropriate  modeling 
of  the  fluid  fluctuations  in  the  near-wall  region  is  crucial  for  accurate  sim¬ 
ulations  of  the  particle  deposition  on  a  surface.  It  is  known  that,  the  k-e 
turbulence  model  is  formulated  implicitly  with  an  isotropic  decomposition  of 
the  turbulence  kinetic  energy  (u'rrns  =  v'rms  =  wrms  =  (2/c/3)1^2).  The  use  of 
(2/c/3)1^2  for  the  velocity  standard  deviation  for  all  components  results  in  an 
overprediction  of  the  normal-to-wall  fluctuating  velocity  component,  which 
is  the  main  component  responsible  for  particle  deposition  on  the  impinge¬ 
ment  wall,  near  the  stagnation  point  and  in  the  wall  jet  region.  As  a  result, 
the  overprediction  of  the  normal-to-wall  fluctuating  velocity  will  cause  an 
overprediction  in  the  deposition  as  was  shown  in  Figures  3.2,  3.3  and  3.4. 

For  this  reason,  some  improvements  are  needed  to  account  for  the  near¬ 
wall  features  in  the  near-wall  region.  Therefore,  to  improve  the  prediction  of 
the  particle  deposition,  a  near-wall  correction  [13]  is  used  in  this  study.  With 
the  near- wall  correction,  the  fluctuating  velocity  in  the  wall  normal  direction 
near  a  wall  can  be  written  as: 

u  =  f  N  u  =  N 

u  JuiyuLLrms 

fu  =  1  -  e~°-02y+ ,  y+ <  80  (3.3) 

where  the  dimensionless  wall  distance  y+  is  calculated  using  the  formula 
y+  =  yu*/u,  and  u*  =  y/rw/p.  Here,  rw  is  the  wall  shear  stress,  and  y  is  the 
normal  distance  to  the  impingement  wall. 

Figures  3.5,  3.6  and  3.7  show  the  deposition  density  curves  as  a  function 
of  r/D  for  L/D  =  2,  4  and  6  obtained  with  turbulent  tracking  and  the  near¬ 
wall  correction.  Table  3.2  shows  the  deposition  efficiency  for  the  three  cases 
obtained  with  turbulent  tracking  and  a  near-wall  correction.  As  evident  in 
Figures  3. 5-3. 7  and  in  Table  3.2,  application  of  turbulent  tracking  with  a 
near-wall  correction  yielded  much  better  results  than  that  without  a  near¬ 
wall  correction  when  compared  to  the  experimental  measurements  [12]. 
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L/D=2 


Figure  3.5:  Deposition  density  for  L/ D  =  2  obtained  from  turbulent  tracking 
and  a  near-wall  correction. 


Figure  3.6:  Deposition  density  for  L/D  =  4  obtained  from  turbulent  tracking 
and  a  near-wall  correction. 


24 


L/D=6 


Figure  3.7:  Deposition  density  for  L/D  =  6  obtained  from  turbulent  tracking 
and  a  near- wall  correction. 


L/D 

Turbulent  tracking 

Near- wall  correction 

CFX  5.7.1 

Exp.  Burwash  (2006) 

2 

56% 

32% 

75% 

16.5% 

4 

57% 

27% 

77.4% 

17% 

6 

41% 

25% 

82% 

17.8% 

Table  3.2:  The  deposition  efficiency  obtained  from  turbulent  tracking  with 
and  a  without  a  near-wall  correction  in  comparison  with  some  experimental 
measurements. 


25 


Chapter  4 

Aerosol  particle  deposition  on 
the  human  form 


4.1  Numerical  parameters 


A  full-size  1.8  meter  tall  human  form  in  the  standing  position  is  shown  in 
Figure  4.1.  The  non-dimensional  computational  domain  size  based  on  the  hu¬ 
man  height  (as  the  reference  length)  are  [—2,  8]  in  ^-direction  (or,  streamwise 
direction),  [—2,2]  in  ^-direction  (or,  spanwise  direction)  and  [—0.5,2]  in  21- 
direction  (or,  vertical  direction).  The  human  form  is  positioned  at  (0,  0,  —0.5) 
in  the  computational  domain  as  illustrated  in  Figure  4.2.  The  human  face 
is  oriented  at  three  orientations  (viz.,  0,  90,  and  180  degrees)  with  respect 
to  the  incident  wind  direction  which  is  assumed  to  be  along  the  +x  direc¬ 
tion.  These  three  orientations  of  the  human  form  with  respect  to  the  incident 
wind  direction  are  shown  in  Figure  4.3  (a),  (b)  and  (c),  which  correspond, 
respectively,  to  0  degree  (human  facing  the  wind),  90  degrees  (human  fac¬ 
ing  direction  90  degree  from  the  incident  wind  direction),  and  180  degrees 
(human  facing  away  from  the  wind). 

The  uniform  inlet  wind  velocity  at  the  x  =  —  2  plane  is  2.5  m  s-1.  The 
Reynolds  number  based  on  the  inlet  wind  velocity  and  height  of  the  human 
form  is 


pUH 

b 


(1.18)  (2.5)  (1.8) 
(1.824  x  lO"5) 


=  2.91  x  105. 


(4.1) 


The  standard  no  slip  boundary  condition  (BC)  is  applied  to  the  walls  (and 
body  surface  of  the  human  form)  and  a  zero  gradient  (Neumann)  BC  com¬ 
bined  with  a  bulk  mass  flow  correction  is  used  at  the  outlet  of  the  computa¬ 
tional  domain.  The  simulation  is  carried  out  in  a  three-dimensional  Cartesian 
framework,  and  the  curved  surfaces  on  the  human  form  that  are  not  aligned 
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Figure  4.1:  A  model  of  the  human  form  used  for  the  simulation  of  aerosol 
particle  deposition. 


with  the  grid  lines  are  approximated  by  stepwise  surfaces.  A  mesh  consisting 
of  100  x  70  x  70  grid  lines  in  the  x-,  y -,  and  ^-directions,  respectively,  is  used 
to  accommodate  all  the  necessary  geometrical  details.  The  grid  arrangement 
adopted  here  is  shown  in  Figure  4.4.  The  grid  lines  are  preferentially  con¬ 
centrated  near  the  human  body  surfaces  (and  other  solid  surfaces)  where  the 
gradients  in  the  flow  properties  are  expected  to  be  greatest,  and  the  spacing 
between  the  grid  lines  is  gently  stretched  with  increasing  distance  from  the 
solid  surfaces. 
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Figure  4.2:  The  non-dimensional  computational  domain  used  for  the  simu¬ 
lation  of  the  aerosol  particle  deposition  on  a  human  form. 
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(c) 

Figure  4.3:  The  human  form  is  directed  at  three  different  orientations  with 
respect  to  the  incident  wind  direction:  (a)  0,  (b)  90,  and  (c)  180  degrees 
corresponding  to  a  human  form  facing  into,  facing  90  degrees  from,  and 
facing  away  from  the  direction  of  the  incoming  wind,  respectively. 
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Figure  4.4:  Generated  three-dimensional  computational  grid. 


4.2  Flow  fields 

The  generated  turbulent  flow  fields  about  the  human  form  for  the  three 
cases,  corresponding  to  the  human  form  oriented  at  0,  90  and  180  degrees 
with  respect  to  the  incoming  wind  direction,  are  summarized  in  Figures  4.5- 
4.10.  For  visualization  purposes,  the  flow  fields  and  streamlines  at  two  cross 
sections  through  the  human  form  are  exhibited:  namely,  in  a  vertical  x-z 
plane  bisecting  the  human  form  along  a  vertical  axis  and  in  a  horizontal  x-y 
plane  bisecting  the  human  form  at  half  height.  From  these  figures,  it  can 
seen  that  complicated  vortices  and  swirling  flow  are  generated  around  and 
behind  the  human  form. 
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Figure  4.5:  Flow  field  (a)  and  streamlines  (b)  in  a  vertical  x-z  plane  of  the 
computational  domain  bisecting  the  human  form  for  the  case  of  a  0  degree 
orientation  relative  to  the  incident  wind  direction. 
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(a) 


(b) 


Figure  4.6:  Flow  field  (a)  and  streamlines  (b)  in  a  horizontal  x-y  plane  of  the 
computational  domain  bisecting  the  human  form  at  half  height  for  the  case 
of  a  0  degree  orientation  relative  to  the  incident  wind  direction. 


32 


(b) 


Figure  4.7:  Flow  field  (a)  and  streamlines  (b)  in  a  vertical  x-z  plane  of  the 
computational  domain  bisecting  the  human  form  for  the  case  of  a  90  degree 
orientation  relative  to  the  incident  wind  direction. 
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(a) 


(b) 


Figure  4.8:  Flow  field  (a)  and  streamlines  (b)  in  a  horizontal  x-y  plane  of  the 
computational  domain  bisecting  the  human  form  at  half  height  for  the  case 
of  a  90  degree  orientation  relative  to  the  incident  wind  direction. 
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Figure  4.9:  Flow  field  (a)  and  streamlines  (b)  in  a  vertical  x-z  plane  of  the 
computational  domain  bisecting  the  human  form  for  the  case  of  a  180  degree 
orientation  relative  to  the  incident  wind  direction. 
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(a) 


(b) 


Figure  4.10:  Flow  field  (a)  and  streamlines  (b)  in  a  horizontal  x-y  plane  of 
the  computational  domain  bisecting  the  human  form  at  half  height  for  the 
case  of  a  180  degree  orientation  relative  to  the  incident  wind  direction. 
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4.3  Aerosol  particle  deposition 

In  this  study,  we  used  thirteen  different  diameters  of  aerosol  particles  in  our 
simulation:  namely,  0.1585,  0.2511,  0.3981,  0.631,  1.0,  1.585,  2.511,  3.981, 
6.31,  10.0,  15.85,  25.118  and  39.812  //.  The  particle  density  is  1008  kg 
m-3.  For  each  particle  size,  one  million  particles  located  with  a  uniform 
random  distribution  on  the  inlet  plane  is  released  (at  x  =  —2)  into  the  non- 
dimensionalized  computational  domain.  The  mesh  representing  the  body 
surface  of  the  human  form  consists  of  178,716  triangular  cell  which  tessellate 
the  entire  body  surface.  The  number  of  particles  that  impact  a  triangular 
cell  of  this  surface  mesh  are  counted  in  each  simulation.  The  deposition  den¬ 
sity  on  each  triangular  cell  is  defined  as  to  be  the  product  of  the  mass  of 
each  particle  and  the  number  of  particles  that  deposit  in  the  cell  divided  by 
the  area  of  the  triangular  cell  with  units  of  g  m-2.  From  this  information, 
the  deposition  probability  in  each  triangular  cell  is  obtained  by  dividing  the 
deposition  density  in  the  cell  by  the  challenge  particle  mass  per  unit  area  on 
the  inlet  plane  (recall  that  the  areal  particle  mass  density  on  the  inlet  plane 
is  uniform). 

The  isopleths  of  the  deposition  probability  for  each  particle  size  and  for 
each  orientation  of  the  human  form  at  that  particle  size  are  exhibited  in 
Figures  4.11-4.49.  In  each  of  these  figures,  the  human  form  is  displayed  in 
four  views:  namely,  front  view,  rear  view,  right  side  view,  and  left  side  view. 
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Figure  4.11:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 


dp  =  39.812  [i m  for  the  human  form  oriented  at  0  degree  with  respect  to  the 


incident  wind  direction. 
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Figure  4.12:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  25.118  /rm  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.13:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  15.85  n m  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.14:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  10.0  n m  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.15:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  6.31  n m  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 


42 


dp=3. 981pm,  0  degree 
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Figure  4.16:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  3.981  pm  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.17:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  2.511  n m  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.18:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  1.585  pm  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.19:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  1.0  /im  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.20:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  0.631  pm  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.21:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  0.3981  n m  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.22:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  0.2511  /rm  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.23:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  0.1585  /rm  for  the  human  form  oriented  at  0  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.24:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  39.812  /mi  for  the  human  form  oriented  at  90  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.25:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  25.118  /im  for  the  human  form  oriented  at  90  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.26:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  15.85  /rm  for  the  human  form  oriented  at  90  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.27:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  10.0  /rm  for  the  human  form  oriented  at  90  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.28:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  6.31  /rm  for  the  human  form  oriented  at  90  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.29:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  3.981  /rm  for  the  human  form  oriented  at  90  degree  with  respect  to  the 
incident  wind  direction. 
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dp=2.511pm,  90  degree 
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Figure  4.30:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  2.511  /rm  for  the  human  form  oriented  at  90  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.31:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  1.585  pm  for  the  human  form  oriented  at  90  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.32:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  1.0  n m  for  the  human  form  oriented  at  90  degree  with  respect  to  the 
incident  wind  direction.. 
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Figure  4.33:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  0.631  /rm  for  the  human  form  oriented  at  90  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.34:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  0.3981  /im  for  the  human  form  oriented  at  90  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.35:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  0.2511  /im  for  the  human  form  oriented  at  90  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.36:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  0.1585  /im  for  the  human  form  oriented  at  90  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.37:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  39.812  /rm  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.38:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  25.118  /rm  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.39:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  15.85  /im  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.40:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  10.0  n m  for  the  human  form  oriented  at  180  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.41:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  6.31  n m  for  the  human  form  oriented  at  180  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.42:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  3.981  /iin  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.43:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  2.511  /im  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction.. 
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Figure  4.44:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  1.585  /im  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.45:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  1.0  /im  for  the  human  form  oriented  at  180  degree  with  respect  to  the 
incident  wind  direction. 
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Figure  4.46:  Isopleths  of  the  deposition  probability  for  a  particle  diameter 
of  dp  =  0.631  n m  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.47:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  0.3981  n m  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.48:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  0.2511  /rm  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 
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Figure  4.49:  Isopleths  of  the  deposition  probability  for  a  particle  diameter  of 
dp  =  0.1585  /rm  for  the  human  form  oriented  at  180  degree  with  respect  to 
the  incident  wind  direction. 


76 


To  provide  the  quantitative  information  for  the  future  use  by  various 
researchers,  we  have  divided  the  human  body  into  eleven  parts  as  follows: 

1.  head  (entire), 

2.  neck  (entire), 

3.  upper  chest  (including  shoulders  and  down  to  just  below  the  nipples), 

4.  lower  abdomen  (including  stomach  and  down  to  the  waist), 

5.  arms  (both  left  and  right), 

6.  hands  (left  and  right), 

7.  legs  (left  and  right,  front  and  back), 

8.  feet  (left  and  right), 

9.  back, 

10.  buttocks, 

11.  crotch  area  (frontal  area  below  the  waist). 

The  decomposition  of  the  human  form  into  the  various  body  parts  along 
with  their  assigned  body  part  indices  is  shown  in  Figure  4.50.  The  following 
information  are  provided  in  the  data  files  that  are  produced  for  each  of  the 
39  simulations  cases  (13  different  particle  diameters  each  associated  with  3 
different  orientations  of  the  human  form  with  respect  to  the  incident  wind 
direction) : 

1.  triangular  cell  number  (index), 

2.  area  of  the  triangular  cell  (m2), 

3.  body  index  associated  with  the  triangular  cell, 

4.  deposition  probability  for  the  triangular  cell. 

A  CD  including  the  data  files  for  each  simulation  case  and  a  C  source  code 
to  read  the  data  sets  from  the  files  is  attached.  The  Appendix  also  includes  a 
listing  of  the  C  source  code.  In  the  CD,  there  are  thirteen  directories  named 
in  accordance  to  the  particle  diameter  (i.e.,  0pl585,  0p2511,  ...,  39p812). 
In  each  directory,  there  are  three  data  files  referenced  by  the  orientation 
angle  of  the  human  form  relative  to  the  incoming  wind  direction  together 
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with  the  particle  size.  For  example,  D_000_0P1585.dat  is  the  name  of  a 
data  file  containing  information  on  the  deposition  probability  on  the  human 
form  for  the  case  associated  with  a  particle  diameter  of  dp  =  0.1585  /irn 
with  the  human  form  oriented  at  0  degree  with  respect  to  the  incident  wind 
direction;  similarly,  D_090_6P31.dat  is  the  name  of  a  data  file  containing 
deposition  probability  on  the  human  form  for  the  case  associated  with  a 
particle  diameter  of  dp  =  6.31  //in  with  the  human  form  oriented  at  90 
degrees  with  respect  to  the  incident  wind  direction. 
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Figure  4.50:  The  human  form  decomposed  into  eleven  body  parts  (and  the 
index  associated  with  each  part). 
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Chapter  5 
Conclusions 


The  URANS  approach  in  conjunction  with  EIM  is  used  here  to  perform 
the  numerical  simulations  of  aerosol  particle  deposition  on  a  realistic  hu¬ 
man  form.  Numerical  simulation  of  a  two-phase  turbulent  impinging  jet  flow 
is  used  to  validate  the  proposed  numerical  methodology  for  the  prediction 
of  the  deposition  density  (or,  equivalently,  the  deposition  probability).  The 
predicted  deposition  results  for  this  case  are  compared  with  some  experimen¬ 
tal  measurements.  The  results  show  that  the  standard  EIM  with  turbulent 
tracking  overpredicts  the  deposition  efficiency.  It  is  demonstrated  that  sig¬ 
nificantly  improved  predictions  for  the  deposition  efficiency  can  be  obtained 
if  a  near- wall  correction  is  included  in  the  EIM. 

After  validation  of  the  numerical  methodology,  detailed  numerical  simu¬ 
lations  of  aerosol  particle  deposition  on  a  human  form  are  conducted.  The 
aerosol  particle  sizes  used  in  these  simulations  cover  the  range  from  about 
0.1  /i m  to  about  40  /mi.  For  each  particle  size,  three  different  orientations  of 
the  human  with  respect  to  the  incident  wind  direction  have  been  simulated 
(namely,  a  human  form  oriented  at  0,  90,  and  180  degrees  with  respect  to 
the  direction  of  the  incoming  wind).  The  detailed  isopleths  of  the  deposi¬ 
tion  probability  on  the  surface  of  the  human  form  have  been  plotted  and  the 
associated  data  sets  have  been  archived  for  future  application.  To  provide 
the  quantitative  information  for  the  future  applications  of  these  data  sets  by 
various  researchers,  the  human  form  has  been  decomposed  into  eleven  body 
parts.  The  association  of  a  body  part  index  with  each  triangular  cell  used  in 
the  tessellation  of  the  surface  of  the  human  form  has  also  been  provided  in 
the  data  files  of  the  deposition  probability. 
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Appendix 


C  source  code  to  read  the  deposition  results  on  the  human: 

^include  <stdlib.h> 

^include  <stdio.h> 

^include  <math.h> 
void  main() 

{ 

/* 

VARIABLES: 

ncells:  the  total  number  of  cells; 
celLarea:  the  area  of  each  cell; 
celLindex: 

1.  head  (entire) 

2.  neck  (entire) 

3.  upper  chest  (including  shoulders  and  down  to  just  below  the  nipples) 

4.  lower  abdomen  (including  stomach  and  down  to  the  waist) 

5.  arms  (both  left  and  right) 

6.  hands  (left  and  right) 

7.  legs  (left  and  right,  front  and  back) 

8.  feet  (left  and  right) 

9.  back 

10.  buttocks 

11.  crotch  area  (frontal  area  below  the  waist) 
depositiomprobability:  deposition  probability  in  each  cell 

7 

double  *cell_area,  *deposition_probability ; 
int  i, id, ncells, *cell_index; 

FILE  *fp; 

fp=fopen(”  human_results.dat” r” ) ; 
fscanf(fp,”  %d”  ,&ncells) ; 

cell_area=  (double*  )malloc(ncells*sizeof  (double)); 
cell_index=(int*)malloc(ncells*sizeof(int)); 
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deposition.probability  =  (double* )  malloc  (ncells*sizeof  (double) ) ; 
for(i=0;  i<ncells;  i++) 

{ 

fscanf(fp,”  %d,%le,%d,%le”  ,&id,&cell_area[i]  ,&cell_index[i]  ,&deposition_probability[i 

} 

fclose(fp); 
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